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We describe an architecture based on a processing 'core' where multiple qubits interact perpet- 
ually, and a separate 'store' where qubits exist in isolation. Computation consists of single qubit 
operations, swaps between the store and the core, and free evolution of the core. This enables com- 
putation using physical systems where the entangling interactions are 'always on'. Alternatively, 
for switchable systems our model constitutes a prescription for optimizing many-qubit gates. We 
discuss implementations of the quantum Fourier transform, Hamiltonian simulation, and quantum 
error correction. 

PACS numbers: 03.67.Lx 



Typically, schemes for solid state quantum computing 
involve an array of qubits with some form of direct phys- 
ical interaction coupling nearby elements [1] . In order to 
implement a specific algorithm, these schemes require the 
experimentalist to dynamically control the magnitude of 
each qubit-qubit interaction - effectively to be able to 
switch it 'on' and 'off'. A common idea for achieving this 
is to somehow dynamically manipulate the wavefunction 
overlap between a pair of neighboring qubits, while other 
nearby qubits are decoupled. This appears feasible, but 
highly challenging. Moreover, even if a switching mecha- 
nism can be implemented, frequent switching is likely to 
increase the rate of dechoerence. A deeper objection is 
that, by having the majority of a system's interactions 
'off' at a given moment, we are failing to maximally ex- 
ploit its computational potential. 

Recently ideas have emerged [2, 3] for computation in 
systems where the interaction remains always on. How- 
ever, these proposals find ways to effectively pacify an 
interaction, and therefore one can make the same ob- 
jection that they are not exploiting the full entangling 
power of the device. One class of system that does make 
full use of a set of permanent interactions is the mirror- 
inversion chain [4, 5, 6, 7, 8, 9]. A chain of spins, with 
suitably engineered coupling strengths, has the property 
that a qubit placed on one end will later manifest at the 
other - even though at intervening times it is distributed 
over the chain. When more than one qubit is placed on 
the chain, each will manifest at the complimentary site - 
but typically the qubits will have aquired an entangling 
phase. It has been observed [7] that this phase could 
in principle be employed to create certain classes of en- 
tangled state, graph states, which are the resource for 
one-way computation. 

In this letter, we demonstrate the potential of such 
engineered spin chains to directly implement arbitrary 
controlled multi-qubit gates. The chain then acts as the 
computation core of our computer (see Figure 1) - we 
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FIG. 1: Architecture of the processor core model. The core 
is an engineered spin- 1/2 chain with always on interactions. 
The storage bank consists of isolated sites where qubits can be 
swapped to and from the corresponding sites in the processing 
core. Controlled multi-target gates are constructed by the free 
evolution of the spin chain. 



need only supplement its free evolution with swap op- 
erations and single qubit manipulations. Note that this 
model is profoundly distinct from schemes involving a 
single qubit bus, e.g. the original ion trap schemes, since 
there the common mode represents only one qubit of in- 
formation. We show that a controlled multi-qubit gate 
can be constructed with exactly four free evolutions of 
the spin chain, independent of the number of spins in- 
volved. The controlling qubit can be any member of the 
spin chain, and the conditional unitary operations ap- 
plied to the target qubits can be of any type. Such a 
gate can significantly reduce the number of elementary 
operations for quantum algorithms involving many non- 
local two-qubit operations. 

We start with a finite chain of N spin-i particles con- 
ned within local potentials and interaction with their 
nearest neighbors. The Hamiltonian considered is: 

JV-l JV 
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where the coupling constants {luj , Xj } are real and in 
general distinct. We will adopt the convention that \0j) 
refers to be the spin-down state ||) (||)) at the 

site j. 
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Next, we require that the core should possess the so- 
called "mirror symmetry" [5] , which implies the coupling 
constants satisfy the relations u>j = lun-j and Xj = Xj 
where j = N—j+1 denotes the mirror-conjugate site of j. 
Let \s) = \s\S2 ■ ■ ■ sat), Sj = {0, 1}, be a particular spin 
configuration. Mirror inversion is said to occur when 
the state s) is driven by the evolution operator U (t) = 
e~ lHr , for some fixed time period r, to the inverted state 
s) = \sn ■ ■ ■ S2S1) (up to a phase factor). Note that the 
term mirror inversion refers to the inversion of classical 
states s) in which all the sites have definite spin values. 
For a quantum state, being a superposition of the basis 
states in general, some internal phases will be acquired. 
However, it is exactly these entangling phases which we 
will exploit for gate construction. 

It is shown [8] that the sufficient and necessary condi- 
tion for mirror inversion in mirror symmetrical chains is 
determined by the eigenvalue spectrum of Hs (h = 1), 
the single excitation subspace of H. Then, 

e~ lEkT = {-if e'^ , (2) 

where Ek, k = 0, 1, 2, ■ • ■ , N — 1, is the (k + l)-th eigen- 
value of Hs and 4>n is some global phase independent of 
k. Due to the mirror symmetry, the coupling constants 
{ujj , Xj } can be determined by the eigenvalue spectrum. 
It is therefore an inverse eigenvalue problem [10]. Re- 
cently, many spectra [4, 5, 9] satisfying the condition in 
Eq.(2) have been proposed. However, to keep our model 
general, we will continue our discussion without reference 
to any specific type of spectrum. 

To construct multi-qubit gates, we need to know the 
matrix elements of the evolution operator U(r) in the 
\s) basis. Let U = U (t). By mapping our picture of 
localized spins to that of spinless fermions [5, 8], one can 
show that 

U |s) = e- nt ^ N (_i)("- m )/ 2 \s) , (3) 

where n is the number of spin-up states in s), and 
?n=0 (1) if n is even (odd). The factor (—i)^ - ™)/ 2 
could be understood intuitively as follows: if n is even 
(i.e. m=0), then the operation of mirror inversion (re- 
ordering the state) is equivalent to swapping n/2 pairs of 
fermions. Similarly for odd n, except the factor should 
be the same as that of n— 1 fermions. The phase factor 
e -m<t> N - IS a coriS equence of Eq.(2). In fact, the phase 4>n 
can be set to zero if an appropriate spectrum of Hs is 
chosen. In this case, it has been demonstrated [7] that the 
operator U (r) alone can generate a fully connected graph 
state. However, in constructing multi-qubit gates, spu- 
rious correlations among qubits in the graph state have 
to be eliminated. This can be achieved with the help of 
an ancilla qubit within the storage array, initialized to be 
|0) a . Let S x represent the swap operation between the 
spin at the site x (compliment of x) and the ancilla. We 
apply S x to the state in Eq.(3) and allow the engineered 



chain to evolve once more, i.e. applying Z x = US X U [11] 
to the initial state |0) Q <E> \s). Then from (3) the final 
state is 

e - { 2n-s^ N g |fiiS2 ...^ ....... y . 

(4) 

which is the same for n being cither odd or even. Here we 
have only assumed a swap operation performed between 
the site x and the ancilla. Therefore, the qubit staying 
at the ancilla spin cannot be transferred back to the spin 
chain at this stage. However, as we shall see (cf. Eq.(7)), 
a more general multi-qubit gate can be constructed based 
on Z x and all of the qubits can reside in their original 
locations at the end of the operation. 

The next step is to interpret the result (4) in terms 
of the quantum circuit model. The phase factor 
e -{2n~s x )uj> N can k e regarded as a result of N local phase 
gates Rj (-20jv), where R,- (cp) = \0 3 ) (0,1 + {l 3 \, 
acting on all qubits, and one extra phase gate, Rj(0at) 
acting on the spin at site x along. On the other hand, 
the factor (—i) 5 ^" -1 ) can be considered as due to the 
application of controlled-tr z to all qubits, except the spin 
at site x which is encoded with the controlling qubit. 
Suppose we now apply local operations to get rid all of 
the controlled phase gates Rj (or simply choose an eigen- 
value spectrum such that 4>n = 0), effectively we have 
constructed a controlled multi-target gate, which requires 
two free evolutions of the engineered Hamiltonian for any 
N. Note that for this multi-qubit gate generated by Z x , 
the o~ z gate has to be applied to all qubits, controlled 
by a single qubit at site x. However, the o~ z gate can be 
converted into controlled- NA, [12], where 

w _ ( sm6j e lVj cos9j\ , , 

' { > ■■' m>«, -sm9j J ' W 

through local operations Aj provided that the relations 
Aj a z Aj = Vj and AjA^ = \j, where \j is the identity 
operator, arc satisfied. For example, the case 6j = (fj=0 
corresponds to a a x gate. Operationally, we denote 
the construction of this controlled multi-target gate by 
V x = AZ X A\ where A = J]jli A i ( and similarly for 
A^). An immediate application of this gate is that, if 
we initialize the controlling qubit to be |0) + |1) and 
the rest |000...0), it can efficiently generate a cat state 
|000...0) + 1 1 1 1 ... 1) , which is interesting for various appli- 
cations including single qubit measurement and encoding 
error correcting codes, such as the Shor's code. 

The controlled operations V 'j are not yet completely 
general: for example, the phase gate Rj((p) and the iden- 
tity operator \j are excluded. We can construct a more 
general controlled multi-target gate, which applies arbi- 
trary unitary operations \Nj on the qubits. To proceed, 
consider applying U to the state \s1S2 • ■ • X • ■ • Sjv), which 
is assumed to contain n spin-up states. The phase fac- 
tors generated are exactly the same as that in Eq.(3). 
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FIG. 2: The circuit diagram of a controlled multi-target 
gate W 1 = AZ 1 BZ 1 C constructed by the free evolution of 
an engineered spin chain 'core'. The parts inside the dashed 
boxes are Z which involves two free evolutions of the spin 
chain and one ancilla in the 'store'. The local operators Aj, 



order to increase efficiency. The primary gain in efficiency 
will of course be a reduction in the number of switching 
events (and consequently, a potential reduction in the 
decohcrence rate) - but remarkably there can also be an 
absolute speed-up by a fixed factor, as we presently dis- 
cuss. 

One of the immediate applications of the controlled 
multi-target gate W 2 is the operation of quantum Fourier 
transform (QFT), which is a key ingredient in many 
quantum algorithms such as the Shor's algorithm. In 
the above notations, the standard QFT circuit can be 
constructed by applying the multi-qubit gates and the 
Hadamard gates alternatively, 



Bj and C, satisfy A^B^C, = lj 



and AjZBjZCj = Wj, where 



the unitary operations Wj, up to a phase, are defined in (6). 



Now we can apply US X to the resulting state; S x re- 
turns \s x ) from the ancilla to the core. The final state is 
e -(2n+s x )i<j) N (^_iy* n w hich can also be considered as 
an controlled multi-target gate analogous to the one gen- 
erated by Z x \s). For simplicity, we again assume 4>n = 0. 
To construct a more general unitary matrix 



QFT = HjvW jv_1 Hjv- 



• •W 2 H 2 W 1 Hi 



(8) 



e *(oy-/3i-«i) CQS7j „ e *(" J -ft+5 J ) sin7j . 
j-\ e i{ aj + f 3 j -S J ) sinlj e i(°y+ft+«:f)coS7j 



(6) 



where Wj = Rj (tt/2 : '~ x ) for j > x and Wj = lj other- 
wise. Here each joint operation Wj costs exactly four free 
evolutions, including two swaps. The QFT circuit depth 
is therefore O(N), as is the total number of switching 
events. For the switched model, the circuit depth is also 
O(N), while the absolute number of switching events is 
0{N 2 ). 

The second application is the simulation of the evo- 
lution of an "artificial" Hamiltonian Ha formally repre- 
senting a joint interaction between r spin-1/2 particles, 



including the identity operator, one can always choose a 



(9) 



set of local operations Aj, Bj and Cj such that the re- which is locally equivalent to the class of the Hamiltonian 



lj and e ia A.ZB.ZC, = Wj [12] are 



lations AjBjCj 
satisfied. Our result is then that the controlled-Wj gate 
\N X , with controlling qubit at site x, can be constructed 
by the following sequence of operations: W x = AZ' T BZ X C, 
where A, B and C are the tensor products of local oper- 



of the form er™ 1 



a,, 



or <7j. Although it is unlikely to find a group of spin- 
1 /2 particles interacting naturally under the Hamiltonian 
Ha, some higher dimensional systems can be mapped by 
those two-level systems. Moreover, the form of Ha can be 
considered as basic building blocks for simulating more 
complex Hamiltonians through local operations and the 
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Bj and Cj 



respectively. Equivalently, 
W x = |0 B )<0 !e |<g)l + |l a )(l,|® JJWj 



(7) 



approximation: 



P i(A+B)At 
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The circuit diagram for an example W 1 is shown in Fig- 
ure 2. For any N, the cost of generating W x includes 
four free evolutions of the engineered chain, two swap 
operations with the same ancilla and local operations. 

We now describe the application of this model to the 
key problems of quantum Fourier transform, Hamilto- 
nian simulation and quantum error correction. When we 
wish to make a statement about the efficiency of our pro- 
cessor core model, we will compare it with a notional 
fully-switched system having a Hamiltonian similar to 
(1), except that the interactions u>% can be independently 
switched on and off. We assume the fully-switched sys- 
tem couples only pairs of qubits simultaneously, although 
it may do so in parallel - i.e. uJiUJi+i = at all times. 
A real fully-switched system could, presumably, activate 
several adjacent interactions: the results we describe here 
can be seen as a prescription for doing precisely that, in 



0(At 2 ) for 
short time-interval At. 

Consider an engineered core of N + 1 spins, with 4>n = 
0, initialized as |0s) = |0siS2...sjv)- The quantum circuit 
for simulating the evolution operator Ua (At) = e ~ lHAAt 
for (9) can be constructed by the following sequence of 
operations [12]: 



U A (At) = H W H To (At) H W H 



(10) 



cr? (or Wj 



where To (At) = exp (— ia^ At) and Wj 
if the qubit at site j is not involved). The basic idea of 
this construction is to store the parity (i.e. m — {0, 1}) 
of the sites j = 1, 2, 3, . . . , N to site 0. The phase gener- 
ated by the local operation To, depending on the parity, 
is exactly the one required for Ha- From (7), it is ap- 
parent that the series of non-local operations W° can be 
achieved by four free evolutions of the engineered chain 
of N > r spins, with the aid of an ancilla and local op- 
erations. This scheme offers a flexibility of generating 
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interactions involving various number of spins using the 
same spin chain. The costs of generating each type of in- 
teractions are fixed (eight free evolutions). Alternatively, 
if one just needs to generate interactions with fixed num- 
ber of qubits, i.e. r = N, the cost can be reduced to two 
free evolutions and no ancilla is needed. The sequence of 
operations in this case is 



Ua (Af) = H UHgTg (At) H UH 



(11) 



The crucial observation for obtaining (11) is that 
H UH |0s) = (_i)("~™)/ 2 \sm). Thus, the desired phase 
can be obtained by applying the local operator Tq at 
site 0. 

The processor core model will also be advantageous 
in running quantum algorithms in a fault tolerant fash- 
ion with concatenated code-words. For example, in the 
Steane code, six gates of the class \N X are required for 
error syndrome measurement (Fig. 10. 16 of Ref.[12]). In 
our approach, each level of concatenation just multiplies 
the number of targets in each Z x by 7 [13] but does not in- 
crease the number of applications of such gates. However, 
the number of elementary switching operations required 
by a fully-switched system to realize \N X increases 7 fold 
with each level of concatenation. 

The discussions above highlight potential gains in 
terms of simplicity: the circuit depth, or total number 
of switching operations. It is also interesting to ask, can 
the total time required be reduced by applying the pro- 
cessor core model? We can quickly conclude that any 
speed-up must be bounded, since the fundamental oper- 
ation U (eqn. 3) can be simulated on a fully-switched 
array in time O(N) [14], while U also takes time O(N) 
to evolve on our processor-core (given a fixed maximum 
interaction strength). Interestingly, there can be speed- 
ups within this constraint, i.e. by bounded factors. To 
make a definite statement we specialize to a core with a 
linear spectrum, i.e. = Ek—Ek+i being constant (e.g. 
[4, 5, 9] ) since this is the time-optimal choice for a given 
spectral range [15]. Let us the compare the time required 
for a simple state transfer, i.e. 1 100. ..0) — ► 1 0. . .001) , 
from one end of the chain to the other. If we now as- 
sume for simplicity that that the maximum interaction 
strength u max ~Y/4 (see Ref. [4]) scales as N, then the 
time required for the processor-core evolution is simply 
ir for all N. On the fully-switched system, the time 
required for each swap is -^j-, thus the total time re- 
quired is T(N) = E^Ti 1 ^"- (Recall that our fully- 
switched system can couple only pairs of qubits simul- 
taneously, thus each state transfer must complete be- 
fore the next is initiated). One can easily show that 
T{N) > (N-l)ir/2u max ~ 2tt (JV-1) /N. Thus this 
always- on processing core can be superior by a factor 
of 2 for large N. 

Finally, we remark that the periods of free processor 
core evolution can be relatively robust versus timing er- 



rors t — > r + St in subsequent swaps to the store. Con- 
sider the most general case, we start with the state \%p) = 
J2j a j where JY |aj| =1. If the evolution time is 
taken perfectly, we expect the final state to be \ip) = 
JY otje l ^i \sj), where <f>j represents the overall phase in 
Eq. (3) for the spin configuration \sj). If not, we have 



(?P\U(T+St) |V> 
1 + iASt + BSt 2 



U{St)\s k ) 



0(St 3 ). 

The error, defined as e = 1 - 
fore just second order in St. 



Here both A and B are real. 

• | (if}\ U (r+Si) |^)| 2 ,is there- 



in conclusion, we have demonstrated how to construct 
controlled multi-target gates through the natural evolu- 
tion of a processor core where interactions are always 
on. This model allows computation with physical sys- 
tems where the entangling interactions are not switch- 
able. Alternatively, in switchable systems our protocol 
can play an important role in simplifying multi-qubit 
operations. We demonstrated this by showing that the 
fully-switched model is fundamentally more complex for 
certain important algorithmic tasks. For various opera- 
tions spanning many qubits, including long range state 
transfer, the temporal requirement for the processor core 
model is less-than-or-equal-to that of the fully-switched 
model. Therefore, in terms of the reduction of dynamical 
control while maintaining the same speed, many schemes 
presently described in terms of two-qubit gates can be 
enhanced by incorporating the processor core concept. 
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